SYSTEM AND ARCHITECTURE FOR PRIVACY-PRESERVING DATA MINING 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The-present invention -relates to-mining-data -from- Internet -users while-preserving the 

privacy of the users. 

2. Description of the Related Art 

The explosive progress in computer networking, data storage, and processor speed has led 
to the creation of very large data bases that record enormous amounts of transactional 
information, including Web-based transactional information. Data mining techniques can then 
be used to discover valuable, non-obvious information from large databases. 

Not surprisingly, many Web users do not wish to have every detail of every transaction 
recorded. Instead, many Web users prefer to maintain considerable privacy. Accordingly, a Web 
user might choose not to give certain information during a transaction, such as income, age, 
number of children, and so on. 

It happens, however, that data mining of Web user information is not only useful to, e.g., 
marketing companies, but it is also useful in better serving Web users. For instance, data mining 
might reveal that people of a certain age in a certain income bracket might prefer particular types 
of vehicles, and generally not prefer other types. Consequently, by knowing the age and income 
bracket of a particular user, an automobile sales Web page can be presented that lists the likely 
vehicles of choice to the user, before other types of vehicles, thereby making the shopping 
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experience more relevant and efficient for the user. Indeed, with the above in mind it will be 
appreciated that data mining makes possible the filtering of data to weed out unwanted 
information, as well as improving search results with less effort. Nonetheless, data mining used 
to improve Web service to a user requires information that the user might not want to share. 

As recognized herein, the primary task of data mining is the development of models about 
aggYegated data. Accordingly, the present invention understands that it is possible to develop 
accurate models without access to precise information in individual data records. Surveys of Web 
users indicate that the majority of users, while expressing concerns about privacy, would willingly 
divulge useful information about themselves if privacy measures were implemented, thereby 
facilitating the gathering of data and mining of useful information. The present invention has 
carefully considered the above considerations and has addressed the noted problems. 

SUMMARY OF THE INVENTION 

The invention is a general purpose computer programmed according to the inventive steps 
herein to mine data from users of the Internet while preserving their privacy. The invention can 
also be embodied as an article of manufacture - a machine component - that is used by a digital 
processing apparatus and which tangibly embodies a program of instructions that are executable 
by the digital processing apparatus to undertake the present invention. This invention is realized 
in a critical machine component that causes a digital processing apparatus to perform the 
inventive method steps herein. The invention is also a computer-implemented method for 
undertaking the acts disclosed below. 
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Accordingly, a computer-implemented method for obtaining data from at least one user 
computer via the Internet while maintaining the privacy of a user of the computer includes 
perturbing original data associated with the user computer to render perturbed data. The method 
also includes generating at least one data mining model using the perturbed data. 

In a preferred embodiment, perturbed data is generated from plural original data associated 
with respective plural user computers. As intended by the present invention, the original data 
cannot be reconstructed from the respective perturbed data. The data can perturbed using a 
uniform probability distribution or a Gaussian probability distribution. Categorical data is 
perturbed by selectively replacing the data with other values based on a probability. 

In another aspect, a computer system includes a program of instructions that include 
structure to, at a user computer, randomize at least some original values of at least some nimieric 
attributes to render perturbed values. The program also sends the perturbed values to a server 
computer, where the perturbed values are processed to generate at least one classification model. 

In still another aspect, a computer storage device includes computer readable code that is 
readable by a server computer for generating at least one classification model based on original 
data values stored at plural user computers without knowing the original values. The device 
includes logic means for receiving perturbed values from the user computers. In accordance with 
the present invention, the perturbed values represent randomized versions of the original values. 
Logic means then generate a classification model using the perturbed values without using the 
original values. 

In yet another aspect, a computer storage device includes computer readable code readable 
by a user computer for facilitating the generation of at least one classification model based on 
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original data values stored at the user computer without knowing the original values. The device 
includes logic means for generating perturbed values representing randomized versions of the 
original values, and logic means for sending the perturbed values to a server computer for 
generating at least one classification model based thereon. 

The details of the present invention, both as to its structure and operation, can best be 
understood in reference to the accompanying drawings, in which like reference numerals refer 
to like parts, and in which: 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a block diagram of the present system; 

Figure 2 is a schematic diagram of a computer program product; 

Figure 3 is a flow chart of the overall logic; 

Figure 4 is a flow chart of the logic for reconstructing the data distribution of the original 
user data; 

Figure 5 is a flow chart of the logic for generating a decision tree classifier; and 
Figure 6 is a flow chart of the logic for generating a Naive Bayes classifier. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Referring initially to Figure 1, a system is shown, generally designated 10, for mining data 
from plural user computers 12 (only a single user computer 12 shown in Figure 1 for clarity of 
disclosure) such that computer-implemented Web sites 14 can more effectively serve the user 
computers 12 while preserving the privacy of the user computers 12. The user computer 12 
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includes an input device 16, such as a keyboard or mouse, for inputting data to the computer 12, 
as well as an output device 18, such as a monitor, for displaying Web pages that have been 
tailored for the particular user of the computer 12. The Web pages are sent via the Internet from, 
e.g., the Web site 14. 

One or both of the computer 12AVeb site 14 can be a personal computer made by 
International Business Machines Corporation (IBM) of Armonk, N.Y. Other digital processors, 
however, may be used, such as a laptop computer, mainframe computer, palmtop computer, 
personal assistant, or any other suitable processing apparatus. Likewise, other input devices, 
including keypads, trackballs, and voice recognition devices can be used, as can other output 
devices, such as printers, other computers or data storage devices, and computer networks. 

In any case, the processor of the user computer 12 accesses a perturbation module 20 to 
undertake certain of the logic of the present invention, while the Web site 14 accesses a privacy 
module 22 to undertake certain of the present logic. The modules 20, 22 may be executed by 
a processor as a series of computer-executable instructions. The instructions may be contained 
on a data storage device with a computer readable medium, such as a computer diskette 24 shown 
in Figure 2 having a computer usable medium 26 with code elements A-D stored thereon. Or, 
the instructions may be stored on random access memory (RAM) of the computers, on a DASD 
array, or on magnetic tape, conventional hard disk drive, electronic read-only memory, optical 
storage device, or other appropriate data storage device. In an illustrative embodiment of the 
invention, the computer-executable instructions may be lines of JAVA code. 

Indeed, the flow charts herein illustrate the structure of the logic of the present invention 
as embodied in computer program software. Those skilled in the art will appreciate that the flow 
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charts illustrate the structures of computer program code elements including logic circuits on an 
integrated circuit, that function according to this invention. Manifestly, the invention is practiced 
in its essential embodiment by a machine component that renders the program code elements in 
a form that instructs a digital processing apparatus (that is, a computer) to perform a sequence 
of function steps corresponding to those shown. 

Now referring to Figure 2, at block 28, the perturbation module 20 perturbs original data 
that the user of the user computer 12 wishes to remain private. For example, the user's age, 
income, and number of children might be perturbed at block 28. In one preferred embodiment, 
the data is perturbed using randomization. 

For numerical attributes such as age and salary, a perturbed value of x^ + r is returned, 
where r is a random value selected from a distribution. In one embodiment, the distribution is 
uniform, i.e., r has a uniform probability distribution between [-a, +q;] with a mean of 0. In 
another embodiment, the distribution is Gaussian, i.e., r has a normal distribution with a mean 
of 0 and a standard deviation a. In contrast, for categorical attributes such as profession, the 
true value of the attribute is returned with a probability p, with a value chosen at random from 
the other possible values for that attribute being returned with a probability of 1-p. 

Proceeding to block 30, in the preferred implementation the perturbed data is sent to the 
privacy module 22 at the Web site 14 via the Internet. Moving to block 32, the privacy module 
22 builds a data mining model, also referred to herein as a classification model, based on the 
aggregated perturbed data from many users. The details of preferred methods for building the 
models, including reconstructing the distribution of the original data, are set forth further below. 
It is noted here, however, that although the preferred method includes reconstructing the 
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distribution of the original data from the distribution of the perturbed data, the Web site 14 does 
not know and cannot reconstruct original data, i.e., the attribute values of individual records from 
any user computer. 

Once a data mining model is generated, several options are possible. For example, at 
block 34 the model can be sent as a JAVA applet to a user computer 12, which can then run the 
model "at block 3 6" oh its original records to determine a classification in accordance with the 
model. For example, the model might determine, based on the user's age and salary and 
assuming that the Web site is, e.g., the site of a vehicle vendor, that the user is of a classification 
that is inclined to purchase sports utility vehicles. The classification, but not the original data, 
can be returned to the Web site 14, which can then send a Web page that has been customized 
for the user's particular classification to the user computer 12 at block 38 for display of the page 
on, e.g., the monitor 18. Accordingly, the returned Web page might display and list SUVs more 
prominently than other vehicles, for the user's convenience, without compromising the privacy 
embedded in the original data, which is not available to the Web site 14. 

Another option is shown at block 40 in Figure 3. If the user has generated a search 
request, the Web site 14 can return to the user the complete search results, along with a data 
mining model for ranking search results based on classification. The user computer 12 can then 
use the model to process its original data to return a classification, which is then used to rank the 
search results as a convenience for the user. Again, however, the user's original data remains 
unavailable to the Web site 14. 

In the preferred embodiment, the data mining model is generated not from a distribution 
of the perturbed data, but from an estimate of the distribution of the original data that is 
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reconstructed from the distribution of the perturbed data, to improve the accuracy of the model. 
The estimate of the original distribution is referred to herein as the reconstructed distribution. 
Figure 4 shows the presently preferred method for generating the reconstructed distribution. As 
noted further below in reference to Figures 5 and 6, the algorithm shown in Figure 4 can be used 
prior to constructing the classification models or it can be integrated into the model generation 
process. Less preferably, in addition to or in lieu of generating the model, if desired the 
reconstructed data can be used for clustering or simply to gain an insight into the profile of the 
users of the system. 

Commencing at block 42, a default uniform distribution is initially assumed, and at block 
44 an integration cycle counter "j" is set equal to zero. Moving to block 46, the derivative of the 
posterior density function f^j^^ can be determined for each attribute "a" using the following 
equation: 

f/i(a):= (1/n) Z (over i=l to n) of { [fy(Wi-a)f^j(a)]/ j (from -oo to +oo) of [fy(Wi.z)f^j(z)dz]}, where 

fx = derivative of the posterior distribution function for the reconstructed 
distribution, fy = derivative of the posterior distribution function for the 
distribution of the perturbed data, n = number of independent random variables Y^, 
Y2,. with y^ being the realization of Y^, it being understood that "Y" herein 
was referred to as "r" in the discussion of Figure 3, Wj = (Xj + y^), a = attribute 
under test, and z is an integration variable satisfying, if Y is the standard normal, 
Fy(z) = (l/((27r)-')e-^^*^^'l 

More preferably, to speed processing time, instead of determining the derivative of the 
posterior density function f^^^' at block 46, a partitioning of the domain of original data values 
for each attribute into "m" intervals "I" is assumed, and a probability Pr(XGIp) that an original 
data point "X" lies within an interval Ip of the original distribution is found as follows. First, the 
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distance between z and Wj (or between a and W;) is jfpproximated to be the distance between the 
midpoints of the intervals in which they He. Also/ the density function fx(a) is approximated to 
be the average of the density function in the interval in which the attribute "a" lies. 
With this in mind, / 



Pr'(XGIp) = (1/n) Z (over s=l to m) of {N(IJ x [(fY(m(I,)-m(Ip))Pr(X e Ip))] / [Z(over t=l to m) 
of (fY(m(l3)-m(I^)Pr(XG/))], where 

I(x) is the interval in which 'yx" lies, m(Ip) is the midpoint of the interval Ip, and 
f(Ip) is the average value off the density function over the interval Ip, p=l,...m. 

Using the preferred memod of partitioning into intervals, the step at block 46 can be 
undertaken in O(m^) time. \j is noted that a naive implementation of the last of the above 
equations will lead to a processing time of O(m^); however, because the denominator is 
independent of Ip, the resiilts of that computation are reused to achieve O(m^) time. In the 
presently preferred emboaiment, the number "m" of intervals is selected such that there are an 
average of 100 data paints in each interval, with "m" being bound 10<m<100. 

It is next determined at decision diamond 48 whether the stopping criterion for the 
iterative process disclosed above has been met. In one preferred embodiment, the iteration is 
stopped when th^reconstructed distribution is statistically the same as the original distribution 
as indicated by/a goodness of fit test. However, since the true original distribution is not 
known, the observed randomized distribution (of the perturbed data) is compared with the is 
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compared with the resuh of the current estimation for the reconstructed distribution, and when 
the two are statistically the same, the stopping criterion has been met, on the^intuition that if these 
two are close, the current estimation for the reconstructed distribution is ^Iso close to the original 
distribution. / 

When the test at decision diamond 48 is negative, the in(!egration cycle counter "j" is 
incremented at block 50, and the process loops back to block /6. Otherwise, the process ends 
at block 52 by returning the reconstructed distribution. / 

Now referring to Figure 5, the logic for constructing a decision tree classifier using the 
reconstructed distribution is seen. Commencing at blo^K 54, for reach attribute in the set "S" of 
data points, a DO loop is entered. Moving to blocky56, split points for partitioning the data set 
"S" pursuant to growing the data tree are evaluated. Preferably, the split points tested are those 
between intervals, with each candidate split point being tested using the so-called "gini" index 
set forth in Classification and Regression Trees , Breiman et al., Wadsworth, Belmont, 1984. To 
summarize, for a data set S containing "ly classes (which can be predefined by the user, if 
desired) the "gini" index is given by l-2pp? where pj is the relative frequency of class "j" in the 
data set "S". For a split dividing "S" irito subsets SI and S2, the index of the split is given by: 

index = ni/n(gini(Sl)) + n2/n(gini(S2)), where Uj - number of classes in SI and 
Yij = number of classes in/S2. 

The data points are a^ociated with the intervals by sorting the values, and assigning the 
N(Ii) lowest values to tha/iirst interval, the next highest values to the next interval, and so on. 
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The split with the highest gini index is then used at block 58 to partition the data set into two 
subsets, with the lower intervals relative to the split point being in one subset and the higher 
intervals being in the other. 

Proceeding to decision diamond 60, it is determined, for each partition, whether most 
elements in the partition are of the same class. If they are not, the logic proceeds to block 62 
for each heterogenous partition to loop back to block 56 to further split that partition. Otherwise, 
when all partitions consist of elements most of which are of the same class, the logic prunes the 
tree at block 64 to remove dependence on statistical noise or variation that may be particular only 
to the training data, in accordance with decision tree prune principles set forth in, e.g., Minimum 
Description Length disclosed by Mehta et al. in "A Fast Scalable Classifier for Data Mining", 
Proc. of the Fifth Int'l Conf. on Extending Database Technology , Avignon, France (1996). The 
pruned tree is returned as the classifier. Thus, it is to be appreciated that since the preferred 
embodiment uses reconstructed data derived from the perturbed data, in a general sense the 
perturbed data is used to generate the classifier. 

When using the logic of Figure 5, the reconstruction logic of Figure 4 can be invoked in 
one of three places. First, the reconstructed distribution can be generated for each attribute once 
prior to executing the logic of Figure 5 using the complete perturbed data set as the decision tree 
training set. The decision tree of Figure 5 is then induced using the reconstructed distribution. 

Or, for each attribute, the training (perturbed) data can first be split by class, and then 
reconstructed distributions generated separately for each class, with the decision tree of Figure 
5 being induced using the reconstructed distribution. Yet again, the by-class reconstruction need 
not be done once at the beginning, but rather at each node during the decision tree growth phase. 
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i.e., just before block 56 of Figure 5. We have found that the latter two methods very accurately 
track the original data and at the same time maintain the inability to know any particular original 
attribute value with any meaningful accuracy. For instance, in one experiment using a synthetic 
data generator and a training set of 100,000 records, the true age value for any particular original 

5 record could not be known, with 95% confidence, within an interval any smaller than 60 years. 

- Nonetheless, the classifier generated by the above-disclosed reconstruction and decision tree logic 

very accurately resembled a similar classification model generated using the original data as a 
control. Moreover, we found that using a Gaussian randomizer at block 28 of Figure 3 resulted 
in even better privacy than using a uniform distribution randomizer, and that use of a Gaussian 



ID randomizer decreased the requirement for the reconstruction logic of Figure 4, although 

combining a Gaussian randomizer with reconstruction improved accuracy vis-a-vis generating a 

y] 

data mining decision tree model using uncorrected Gaussian-perturbed data, 
d] Figure 6 shows that as an alternative to generating a decision tree classifier, a Naive Bayes 

= classification model can be generated. Commencing at block 66, the classes C: of data are 

a ^ 

ft determined empirically or using a decision tree-like grovrth phase such as the one shown in 



^1 Figure 5, with the leaf nodes of the tree defining the classes. Moving to block 68, the probability 

Pr(ai = Vy I Cj) of the i* attribute "a" of a record having a value v^ belonging to the class is 
determined by determining the ratio of the number of records in the j'^ class, divided by the total 
number of records, using the perturbed data as a training set. 

20 Next, the logic determines proceeds to block 70, wherein the probability Pr(r | Cj) of a 

record "r" given a class Cj is determined to be 11 (i=l to n) of Pr(ai=Vj | Cj), where aj is an 
attribute that has the value Vj. As before, the preferred way to undertake the above calculation 
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is to partition the perturbed (training) data set into "m" intervals "I", and approximate 
Pr(aj=Vi I Cj) with Pr(aj G P/Cj = the number of records whose class is Cj, where the value of 
the attribute a^ is in the i^ interval partition. Mathematically, this is expressed as {N(Cj A aj G 
Pi)}/N(Cj). Also, the class probability Pr(Cj) of a class occurring is determined using the training 
set. 

After the step at block 70, at block 71 the probability Pr(Cj | r) of a record r being in the 
class Cj is determined by combining the probability Pr(r | Cj) with the class probability Pr(Cj). 
In a preferred embodiment, this is done by multiplying the value found at block 68 by the value 
determined at block 70 (i.e., {N(CjAai G Pi)}/N(Cj). The set of these probabilities for the 
various classes identified at block 66 is then returned at block 72 as the Naive Bayes classifier. 

When using the logic of Figure 6, the reconstruction logic of Figure 4 can be invoked in 
one of two places. First, the reconstructed distribution can be generated for each attribute once 
prior to executing the logic of Figure 6 using the complete perturbed data set as the training set. 
The Naive Bayes classifier of Figure 6 is then induced using the reconstructed distribution. Or, 
for each attribute, the training (perturbed) data can first be split by class, and then reconstructed 
distributions generated separately for each class, i.e., after block 66. 

While the particular SYSTEM AND ARCHITECTURE FOR PRIVACY-PRESERVING 
DATA MINING as herein shown and described in detail is fully capable of attaining the above- 
described objects of the invention, it is to be understood that it is the presently preferred 
embodiment of the present invention and is thus representative of the subject matter which is 
broadly contemplated by the present invention, that the scope of the present invention fully 
encompasses other embodiments which may become obvious to those skilled in the art, and that 



IBM Case No. AM9-99-0226 



- 13 - 



• # 

the scope of the present invention is accordingly to be limited by nothing other than the appended 
claims, in which reference to an element in the singular is not intended to mean "one and only 
one" unless explicitly so stated, but rather "one or more". All structural and functional 
equivalents to the elements of the above-described preferred embodiment that are known or later 
come to be knovm to those of ordinary skill in the art are expressly incorporated herein by 

reference~and^^intended to be encompassed by the present claims. Moreover, it is not 

necessary for a device or method to address each and every problem sought to be solved by the 
present invention, for it to be encompassed by the present claims. Furthermore, no element, 
component, or method step in the present disclosure is intended to be dedicated to the public 
regardless of whether the element, component, or method step is explicitly recited in the claims. 
No claim element herein is to be construed under the provisions of 35 U.S.C. §112, sixth 
paragraph, unless the element is expressly recited using the phrase "means for" or, in the case of 
a method claim, the element is recited as a "step" instead of an "act". 
WE CLAIM: 
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